home *** CD-ROM | disk | FTP | other *** search
/ Shareware Grab Bag / Shareware Grab Bag.iso / 090 / ust1a.arc / FCRC16.C < prev    next >
Encoding:
C/C++ Source or Header  |  1985-09-05  |  4.5 KB  |  160 lines
  1. /*+
  2.  * File: fcrc16
  3.  * Description:
  4.  *    Calculates crc, given a file pointer.
  5.  *
  6.  *
  7.  *  Author: Mohsen Banan.
  8.  *
  9.  *  This program is public domain software, no warranty intended or
  10.  *  implied.
  11.  *  General permission to copy or modify, but not for profit,  is
  12.  *  hereby  granted.
  13.  *
  14.  *
  15.  * Functions:
  16.  *
  17.  *
  18.  * Audit Trail:  $Log:    fcrc16.c,v $
  19.  * Revision 1.1  85/08/28  08:38:33  mohsen
  20.  * original posting to the net
  21.  * 
  22.  * 
  23.  *
  24. -*/
  25.  
  26. #ifdef RCS_VER
  27. static char *rcs = "$Header: fcrc16.c,v 1.1 85/08/28 08:38:33 mohsen Exp $";
  28. #endif
  29.  
  30. /* #includes */
  31. #include "mbstd.h"
  32. #include <stdio.h>
  33.  
  34. /* #defines */
  35.  
  36. /* external variables */
  37.  
  38. /* referenced external function declarations */
  39.  
  40. /* internal function declarations */
  41.  
  42. /* global variables */
  43.  
  44. /* static variables */
  45.  
  46.  
  47. /*<
  48.  * Function:crc16_table
  49.  * Description:
  50.  *
  51.  *  * Calculate CRC16 of a file.
  52.  *  * CRC16 polynomial: x**0 + x**2 + x**15 + x**16 (0xA001)
  53.  *  * (CCITT polynomial: x**0 + x**5 + x**12 + x**16 (0x8408))
  54.  *  * 
  55.  *  * D. Hugh Redelmeier  1983 April 15
  56.  *
  57.  *  Mohsen Banan: I don't know how CRCs work yet. 
  58.  *  None of what is to follow is original work.
  59.  *
  60.  *
  61.  * Returns:
  62.  * 
  63. >*/
  64. #if 0
  65. VOID 
  66. crc16_table()
  67. {
  68.     register unsigned s, t;
  69.     register int i;
  70.  
  71.     printf("unsigned short crctab[256] = {");
  72.     for (t=0; ; ) {
  73.         if (t%8 == 0)
  74.             printf("\n");
  75.         s = t;
  76.         for (i=8; i--; )
  77.             s = s>>1 ^ (s&1? 0xA001 : 0);
  78.         printf("\t0x%04x", s);
  79.         if (++t == 256)
  80.             break;
  81.         printf(",");
  82.         }
  83.     printf("};\n");
  84. }
  85. #endif
  86.  
  87. unsigned short crctab[256] = {
  88.     0x0000, 0xc0c1, 0xc181, 0x0140, 0xc301, 0x03c0, 0x0280, 0xc241,
  89.     0xc601, 0x06c0, 0x0780, 0xc741, 0x0500, 0xc5c1, 0xc481, 0x0440,
  90.     0xcc01, 0x0cc0, 0x0d80, 0xcd41, 0x0f00, 0xcfc1, 0xce81, 0x0e40,
  91.     0x0a00, 0xcac1, 0xcb81, 0x0b40, 0xc901, 0x09c0, 0x0880, 0xc841,
  92.     0xd801, 0x18c0, 0x1980, 0xd941, 0x1b00, 0xdbc1, 0xda81, 0x1a40,
  93.     0x1e00, 0xdec1, 0xdf81, 0x1f40, 0xdd01, 0x1dc0, 0x1c80, 0xdc41,
  94.     0x1400, 0xd4c1, 0xd581, 0x1540, 0xd701, 0x17c0, 0x1680, 0xd641,
  95.     0xd201, 0x12c0, 0x1380, 0xd341, 0x1100, 0xd1c1, 0xd081, 0x1040,
  96.     0xf001, 0x30c0, 0x3180, 0xf141, 0x3300, 0xf3c1, 0xf281, 0x3240,
  97.     0x3600, 0xf6c1, 0xf781, 0x3740, 0xf501, 0x35c0, 0x3480, 0xf441,
  98.     0x3c00, 0xfcc1, 0xfd81, 0x3d40, 0xff01, 0x3fc0, 0x3e80, 0xfe41,
  99.     0xfa01, 0x3ac0, 0x3b80, 0xfb41, 0x3900, 0xf9c1, 0xf881, 0x3840,
  100.     0x2800, 0xe8c1, 0xe981, 0x2940, 0xeb01, 0x2bc0, 0x2a80, 0xea41,
  101.     0xee01, 0x2ec0, 0x2f80, 0xef41, 0x2d00, 0xedc1, 0xec81, 0x2c40,
  102.     0xe401, 0x24c0, 0x2580, 0xe541, 0x2700, 0xe7c1, 0xe681, 0x2640,
  103.     0x2200, 0xe2c1, 0xe381, 0x2340, 0xe101, 0x21c0, 0x2080, 0xe041,
  104.     0xa001, 0x60c0, 0x6180, 0xa141, 0x6300, 0xa3c1, 0xa281, 0x6240,
  105.     0x6600, 0xa6c1, 0xa781, 0x6740, 0xa501, 0x65c0, 0x6480, 0xa441,
  106.     0x6c00, 0xacc1, 0xad81, 0x6d40, 0xaf01, 0x6fc0, 0x6e80, 0xae41,
  107.     0xaa01, 0x6ac0, 0x6b80, 0xab41, 0x6900, 0xa9c1, 0xa881, 0x6840,
  108.     0x7800, 0xb8c1, 0xb981, 0x7940, 0xbb01, 0x7bc0, 0x7a80, 0xba41,
  109.     0xbe01, 0x7ec0, 0x7f80, 0xbf41, 0x7d00, 0xbdc1, 0xbc81, 0x7c40,
  110.     0xb401, 0x74c0, 0x7580, 0xb541, 0x7700, 0xb7c1, 0xb681, 0x7640,
  111.     0x7200, 0xb2c1, 0xb381, 0x7340, 0xb101, 0x71c0, 0x7080, 0xb041,
  112.     0x5000, 0x90c1, 0x9181, 0x5140, 0x9301, 0x53c0, 0x5280, 0x9241,
  113.     0x9601, 0x56c0, 0x5780, 0x9741, 0x5500, 0x95c1, 0x9481, 0x5440,
  114.     0x9c01, 0x5cc0, 0x5d80, 0x9d41, 0x5f00, 0x9fc1, 0x9e81, 0x5e40,
  115.     0x5a00, 0x9ac1, 0x9b81, 0x5b40, 0x9901, 0x59c0, 0x5880, 0x9841,
  116.     0x8801, 0x48c0, 0x4980, 0x8941, 0x4b00, 0x8bc1, 0x8a81, 0x4a40,
  117.     0x4e00, 0x8ec1, 0x8f81, 0x4f40, 0x8d01, 0x4dc0, 0x4c80, 0x8c41,
  118.     0x4400, 0x84c1, 0x8581, 0x4540, 0x8701, 0x47c0, 0x4680, 0x8641,
  119.     0x8201, 0x42c0, 0x4380, 0x8341, 0x4100, 0x81c1, 0x8081, 0x4040};
  120.  
  121.  
  122.  
  123. /*<
  124.  * Function:fcrc16
  125.  * Description:
  126.  *  The algorithm comes from 
  127.  *  "Byte-wise CRC Calculations" by Aram Perez,
  128.  *  IEEE Micro June 1983 pp.40-50
  129.  *
  130.  * Returns:
  131.  *  VOID
  132.  * 
  133. >*/
  134.     
  135. PUBLIC USHORT 
  136. fcrc16 (i_fp, prev_crc, p_fsize)
  137. FILE * i_fp;
  138. USHORT prev_crc;
  139. LONG * p_fsize;
  140. {
  141.     USHORT crc, x;
  142.     INT c;
  143.     LONG  fsize;
  144.  
  145.     crc = prev_crc;
  146.     fsize = 0L;
  147.     while ((c=getc(i_fp)) != EOF) {
  148.         /* these are really the same.
  149.          * crc = crctab[(crc^c)&0xff] ^ crc>>8;
  150.          */
  151.         x = (c^crc) & 0xFF;
  152.         crc >>= 8;
  153.         crc ^= crctab[x];
  154.         ++fsize;
  155.     }
  156.     rewind(i_fp);        /* I don't like side effects */
  157.     *p_fsize = fsize;
  158.     return (crc);
  159. }
  160.